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Predefinition of gene expression 
rate wr preoeterminea gene 
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Determination of resultant gene expression 
pattern using causal network 
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Comparison of resultant gene expression 



pattern <-> predefined gene 
expression pattern 
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Sampling algorithm (B,N) 
Input: 

B - Bayesian network 
N - Number of independent samples 
Output 
DB- 



1. Arrange the variable set X according to the condition 
so that parents Paj are arranged before Xj 

2. For s=1 N-230 

3. Fori=1,...,n - 220 

4. Let Xj be the node with the highest sequence number 
in this sample which is not instantiated 

5. If Xj is a root node, select the state with 



-210 



6. Otherwise select the state with probability 
P(state|sampled states of PaJ. 

7, Instantiate XNstate 
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Interventional sampling algorithm (B,E,N) 
Input: 

B - Bayesian network 
E- Set of observations 
N - Number of independent samples 
Output: /330 
D B | E - data set of N independent samples for given E. 

XE - Set of observed variables; 320 
Xq={X\XE} - Set of request varaibles 

1. Arrange Xq according to the condition that 
parents Pas are arranged before Xj 

2. For s=1 N 

3. For i=1 n 

4. Let Xj be the node with the highest sequence number 
in this sample which is not instantiated 

5. If Xj is a root node, select the state with 
probability P(state\E) 

6. Otherwise select the state with probability 
P(state|sampled states ofPajE). 

7. Instantiate Xj=state 
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REPLACEMENT SHEET (RULE 26) 



